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I The problem of time reparameterization is addressed at both the classical and quantum levels 

I in a Bianchi-I universe in which the matter source is a massive Dirac spinor field. We take the 

I scale factors of the metric as the intrinsic time and their conjugate momenta as the extrinsic 
time. A scalar character of the spinor field is identified as a representation of the extrinsic time. 

^ ' The construction of the field equations and quantization of the model is achieved by solving the 

<sJ '< Hamiltonian constraint after time identification has been dealt with. This procedure leads to a 

^ ' true Hamiltonian whose exact solutions for the above choices of time are presented. 
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^ . 1 Introduction 
00 

G^ 

^SJ ' Standard cosmological models based on classical general relativity have no convincing precise answer 

i for the presence of the so-called "Big-Bang" singularity. Any hope of dealing with such singularities 

' would be in vein unless a reliable quantum theory of gravity can be constructed. In the absence of 

■ a full theory of quantum gravity, it would be useful to describe the quantum states of the universe 

within the context of quantum cosmology, introduced in the works of DeWitt [1] and later Misner 



[2]. In this formalism which is based on the canonical quantization procedure, one first freezes a large 
^ ■ number of degrees of freedom and then quantizes the remaining ones. The quantum state of the 

universe is then described by a wave function in the mini-superspace, a function of the 3-geometry of 
the model and matter fields presented in the theory, satisfying the Wheeler-DeWitt (WD) equation. 
In more recent times such works have been the focus of an active area of research with different 
approaches, [3]-[ll], see also [12] for a review. In references [13], canonical quantization is applied to 
many models with different matter fields as the sources of gravity. 

As are well known, quantum cosmology suffers from a number of problems, namely the con- 
struction of the Hilbert space to define a positive definite inner product of the solutions of the WD 
equation, the operator ordering problem and also most importantly, the problem of time. The wave 
function in the WD equation is independent of time, i.e. the universe has a static picture in this 
scenario. This problem was first addressed in [1] by DeWitt himself. However, he argued that the 
problem of time should not be considered as a hinderance in the sense that the theory itself must 
include a suitable well-defined time in terms of its geometry or matter fields. In this scheme time is 
identified with one of the characters of the geometry, usually the scale factors of the geometry and 
is referred to as the intrinsic time, or with the momenta conjugate to the scale factors, or even with 
a scalar character of matter fields coupled to gravity in any specific model, known as the extrinsic 
time. 
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In general, the crucial problem in canonical quantum gravity is the presence of constraints in the 
gravitational field equations. Identification of time with one of the dynamical variables depends on the 
method we use to deal with theses constraints. Different approaches arising from these methods have 
been investigated in detail in [14]. The issue of time in canonical general relativity is also extensively 
discussed in [15]. As has been discussed in [14], time may be identified before or after quantization 
has been done. There are approaches, on the other hand, in which time has no fundamental role. 
For a more modern review of the problem of time and other related problems in quantum cosmology 
see [16]. The details of time identification procedure in terms of various dynamical variables of the 
theory before quantization is done has been investigated in [17] where a Robertson- Walker universe 
filled with a scalar field is quantized. Also in [18] a choice of time in terms of a massless scalar field 
is discussed in a Bianchi-I classical cosmology based on the method developed in [19]. 

One of the common candidates of time in the above works is the matter field present in the theory, 
that is, time is identified with a scalar character of matter. In the case of a scalar field as the source 
of gravity, the scalar field itself can play the role of time as is the case in [17] and [18]. Another 
matter field which has occasionally been studied in the literatures is the massless or massive spinor 
field as the source of gravity. In general, theories studying spinor fields coupled to gravity result in 
Einstein-Dirac systems which arc not easy to solve. The quantized Robertson- Walker or Bianchi-I 
universe filled with a spinor field are studied in [20]- [25]. For a general discussion on the possibility 
that classical homogeneous spinor fields might play the role of matter in cosmology, the reader is 
referred to [26]. 

In this paper we deal with classical and quantum cosmology of a model in which a classical massive 
spinor field is coupled to gravity in a Bianchi type I space-time. What we mean by a classical spinor 
field is a set of four complex-valued space time functions that transform according to the spinor 
representation of the Lorentz group. The existence of such fields is crucial in our work since in spite 
of fact that fermions are described by quantized spinor fields which do not have a classical limit, 
we assume such classical fields exist and use them as matter fields coupled to gravity. A possible 
justification for the existence of classical spinor fields is given in the appendix of reference [26]. To 
identify time, we have adopted the same procedure as in [17], that is, after choosing a time parameter, 
we solve the Hamiltonian constraint equation to obtain a minimal true Hamiltonian. We then use 
this reduced Hamiltonian to construct the classical field equations. Also, to quantize the model we 
use the operator form of the resulting Hamiltonian in the Schrodinger picture. The reason that we 
have used the Schrodinger picture for quantizing the system and not the WD picture is that time 
naturally enters the Schrodinger equation. As is well known in the WD formalism, in quantizing a 
cosmological system, there is no a priori definition of time. This is quite natural and stems from the 
nature of the WD equation and is a reflection of the gauge invariance with respect to the choice of 
coordinates in the classical theory. The intrinsic times which we have chosen in this paper are related 
to the scale factors of the metric, with their conjugate momenta chosen as the extrinsic time. We 
have shown that there is also an identification of the extrinsic time in terms of the scalars constructed 
with spinor fields. The exact solutions of the classical and quantum cosmology corresponding to each 
choice of time are presented. The complicated form of some of these solutions show that the evolution 
of the universe with respect to diflFerent intrinsic or extrinsic times is a non-trivial undertaking. 

2 The classical model 

We start with a space-time metric of the form ^ 

= -N\t)dt^ + a^{t)dx'^ + + c^(t)dz^^ (1) 

which describe a Bianchi type I universe with scale factors a{t), b{t) and c(t) in the x, y and z 
directions respectively, with N(t) being the lapse function. This metric is the simplest anisotropic 

^We work in the units where c = % = IQnG = 1 
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and homogeneous cosmological model which, upon making the scale factors equal, becomes the flat 
Robertson- Walker metric. Such space-times have an Abelian symmetry group of translations with 
Killing vector fields ^ = {dx, dy, dz)- Of course, all the structure constants of such a symmetry group 
are zero. The properties of such a space time is studied in many works, see for example [22] and [27] 
and the references therein. The scalar curvature corresponding to metric (1) is 

2 /a h c ah he ca Na Nh Nc\ 
~l^\^^h^~c^d)^h'c^m~l^~Wb~Wc)' ^' 

where a dot represents differentiation with respect to t. To construct the field equations, let us start 
with the action 

^ ^ j (^grav + Lmatt)V^d'^X, (3) 

where 

Lgrav =n-2A, (4) 

is the Einstein-Hilbert Lagrangian for the gravitational field with cosmological constant A, and Lmatt 
represents the Lagrangian of the matter source which we assume to be a classical massive free spinor 
field. As we mentioned in introduction, for the bulk of this paper we shall consider the spinor fields 
as classical objects, i. e. four complex-valued space time functions and not a set of Grassmanian 
variables. In fact spinors in quantum field theory are four operator- valued fields which act on the 
corresponding Hilbert space and obey the Dirac equation (in fiat space time) 

{-yf'd^ - m)il) = 0. (5) 

However, in what follows, following [26], we shall interpret a classical spinor field, that is, a set of 
four complex- valued functions as 

11)^1 = {s \ ip \ s), (6) 

where | s) is an appropriate physical state. Now, we can see that the expectation value of a spinor in 
a physical state is a complex number and not a Grassmanian number. Also taking the expectation 
value of the Dirac equation (5) yields 

(7''^^ - m)Ai = 0, (7) 

which means that the classical spinor fields also obey the Dirac equation. In what follows, by we 
mean the classical field ipd and omit the subscript cl from now on. For a more extensive discussion 
of the properties of such classical fields see [26] . 

As usual, the Dirac equation describing the dynamics of a spinor field can be obtained from the 
Lagrangian 



■^matt — 2 



where 7^ are the Dirac matrices associated with the space-time metric satisfying the Clifford algebra 
{t'^jT'^} = 2^^'', are spin connections and V{ijj,ip) is a potential describing the interaction of the 
spinor field with itself. In the case of a free spinor field of mass m we have ¥{1/;, ip) = mtptp. The 7^^ 
matrices are related to the flat Dirac matrices, 7", through the tetrads as follows 

r = e^r, 7m = <7a- (9) 

For the metric (1) the tetrads can be easily obtained from their definition, that is g^iy = e^^e^rjab, 
leading to 

e^l = diag(iV, a, b, c), = diag(l/A^, 1/a, 1/6, 1/c). (10) 
Also, the spin connections satisfy the relation 

r^ = ^5.A(Va + r^M<hV- (11) 
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Thus, for the line element (1), use of (9) and (11) yields 



o„i 



ro = o, ri = -— 7^ 



2 -n - u 



0,,3 



2N 



2N 



(12) 



Here 7^^ and 7* are the Dirac matrices in Minkowski space and we have adopted the following repre- 
sentation [28] 

n f -i 0\ f a' \ 



7 



i 



7 



a' 



(13) 



The final remark about Lagrangian (8) is that consistency of Einstein field equations with a spinor 
field as the matter source in the background metric (1) requires the spinor field ^ to be dependent 
on t only, that is i/j = ip{t) [24]. 

The preliminary set-up for writing the action is now complete. Substituting (2), (4) and (8) 
into (3) and integrating over the spatial dimensions, we are led to an effective Lagrangian in the 
mini-superspace {N, a, b, c, i/^, ijj}^ 



jC= — (abc + abc + abc) + ANabc + -Nabc 
N \ J 2 



(V'7°'0-'07V)-2F(V','0) 



(14) 



Variation of Lagrangian (14) with respect to ip, ^, a, b, c and A*" yields the equations of motion of 
the spinor and the gravitational fields as 



1 f d b c 



dV 



+ 7 + - U - ^^^7 = 0, 



2 \ a b c 



b c be N b c\ , 1 
b c be N \ b c 2 



c a ac N f c a 

- + - + TT - + - 



c a ac N \c a 



d b db N I d b 



AiV2 = -Af2 



+ T + 



+ 



a b ab N \ a b 



-dV _^ dV^ 
d-i/j dip 

dip dip 
jdV dV ' 



V{ip,ip), 
-V{iP,iP), 
■V{iP,iP), 



db dc be . Ar2T^/ 7 ,x 

— + — + — - AN^ = -N^V{ip,ip). 
ab ac be 



(15) 

(16) 
(17) 
(18) 
(19) 
(20) 



As a double check, one may obtain the above field equations from the Dirac and Einstein equations, 
given as 

(21) 



(22) 



(23) 



^Although, it is not a priori evident that the substitution of the ansatz for the metric and the matter fields into the 
action leads to the correct equations of motion, here, as we can see from the equations, this is the case. This procedure 
is not correct when class B Bianchi models are involved. 
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For a homogeneous spinor field = equations (21) and (22) are equivalent to (15) and (16) 

respectively. On the other hand, the non-vanishing components of the Einstein tensor for metric (1) 
are 

^ ab ac be 

Goo = — + — + -, 24 

ab ac be 

_ a% a^c a^bN a^cN a%c 

with the cyclic permutations a — 6 — c on Gn giving the expressions for G22 and G33. Also the 
components of the energy-momentum tensor for the spinor field as the matter source can be obtained 
from the standard definition 

yielding 

Too = -N^V{^,^), rn = ia2(^V^|I + |^v) T,,- = Tq^ = 0. (26) 

Again the above mentioned permutations on Tn lead to T22 and T^^. Substitution of these results 
into Einstein equations (23) yields the same equations as (17)-(20). 

Let us now construct the Hamiltonian for our model. The momenta conjugate to the dynamical 
variables are 

dC 1 , dC 1 , . ., dC 1 , ; 



and 



dC 1 , T _ _ 1 J, 



p^ = —j = -abcipj , pT = — = --abcy ip. (28) 
dtp 2 2 



The primary constraints are given by 

d/L 1 1 

PN = —y- =0, 7^4,= Pi, - -abcijj-f^ = 0, TTj; = pj; -abcj^ip = 0. (29) 

In terms of the conjugate momenta the Hamiltonian is given by 

H = paCL + pbb + PcC + p^ip + ipp^ + pnN - C, (30) 

leading to 

N f a 2 , b 2 , c 2\ , N ( Pa,pb , PbPc , PaPc 



Because of the existence of constraints (29), the Lagrangian of the system is singular and the total 
Hamiltonian can be constructed by adding to H the primary constraints multiplied by arbitrary 
functions of time A(i), A^(t) and A^(t) 

N f a 2 b 2 , ■7\ , PaPb , PbPc , PaPc 

Ht = --r T-Pa + —Pb + -rPc j+Tr[ + + —r- 

4 \bc ac ab J 2 \ c a b 

- ANabc + NabcV{4>, ip) + XpN + 7rv,A^ + A^tt^, (32) 

where the subscript T stands for total. The requirement that the primary constraints should hold 
during the evolution of the system means that 

p'n = {PN, Ht} ~ 0, (33) 
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and 

tt'v, = {vr^, Ht} ~ 0, 7f^ = Ivr^, ifr} ~ 0. 
The condition (33) leads to the secondary constraint 



(34) 



1 / a 



4 \bc 



ac 



ah 



1 fPaPh . PbP. 



+ 



+ 



PaPc 



- Aabc + abcV {i),i>)=0, (35) 



while relations (34) only fix the functions and 

A^ - 



(36) 



The study of the algebra of constraints (29) together with the Poisson brackets of the secondary con- 
straint with other conjugate variables show that there are no additional constraints in this dynamical 
system [29]. 

Although the field equations obtained above can be solved after a suitable form for the potential 
V{'ip,ip) has been chosen [30, 31], the Hamiltonian constraint (35) does not have the desired form for 
lh(^ construction of the quantized model or for the discussion of the problem of time in this model. 
Thus, to transform Lagrangian (14) to a more manageable form, consider the following change of 
variables 



u+v— 



h = e"" " " , c = e 
In terms of these new variables, Lagrangian (14) takes the form 



u—2v 



AT V / 2 



The momenta conjugate to u, v and w are 



dC 6 



Pv 



-ue 



3« 



dC 



Pv 



dii N dv 
giving rise to the following Hamiltonian 

N 
12' 



N 



dC 
dw 



N 



-we 



3u 



with the corresponding Hamiltonian constraint becoming 



n 



e^" = 0. 



(37) 



(38) 



(39) 



(40) 



(41) 



Now, variation of Lagrangian (38) with respect to its dynamical variables yields the following field 
equations 

y,+ ^^^ + iV/|I = 0, 



■- 3 - BV n ^ 



— .. 3.9 „ N . 3 , . 9 . 9^ 



2 
iV 



KN - -N 
2 



l(Vi7V-V^7»-2mV') 



0, 



N 



-ve 



3u 



3« 



= 0, 
= 0, 



(42) 
(43) 

(44) 
(45) 
(46) 
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3 (.2 



w^) + - A) =0. 



(47) 



Up to this point the cosmological model, in view of the concerning issue of time, has been rather 
general and of course undcr-dctcrmined. Before trying to solve these equations we must decide on a 
choice of time in the theory. The general solutions of the system (42)- (47) for a free spinor field can 
be written as [25] 



u{t) = In 



m 



m 
^2A 



+ 



+ 



4A2 



+ 



-3A 



cos 



V- 



3A / N{t')dt' 
'to 



+ 



C2 



4A2 -3A 



cos 



3A f N{t')dt' 



1/3 



v{t) = I In 



w{t) = ^ In 



{B - ^)tan (^v^=3A//^ N{t')dt') + C/^T^ 



{B - A) tan (i^/^j/^ N{t')dt'^ - Cj^- 
{B - A) tan (i^/^//^^ N{t')dt') + C/^f- 



3A 



3 A 



{B - A) tan (^^V^SA Jl N{t')dt') - C/V^3A 



(48) 



(49) 



(50) 



(51) 



where A = 3^ and B = y ^ + 3^ with ci and C2 being two integrating constants such that 
= 9(cf + C2). To obtain the above solutions we have assume that A < 0. The corresponding 
solutions for A > can easily be obtained by the replacement of the cosine function in (49) with its 
hyperbolic counterpart. The under-determinacy problem at the classical level may be removed by 
using the gauge freedom via fixing the gauge. For example, we can work in the gauge iV = 1 which 
usually is chosen in classical cosmological models and is called the cosmic time gauge. It is also worth 
noting that with the form (1) as the space-time metric, we have done the first step in gauge fixing, 
namely = 0, where is the shift vector. For other candidates in gauge fixing see [14] and [17]. In 
any case, after fixing the gauge, the lapse function N(t) and shift vectors N^(t) are eliminated from 
the field equations, rendering them as the Hamiltonian and momentum constraints. Elimination of 
the lapse function from equations (48)-(51) gives rise to the intrinsic dynamics of the system as a 
relationship between the 3-geometry and matter field which is independent of the choice of time in a 
particular gauge. The result is 

= (52) 



or 



tt(V'V') = In 



1\ 1/3 



'ip'ip 



v{iljtp) 



£1 
C 



In 



B-A 

l+Aipijj 



^ [{B - A)^i; - 1] [{B + A)i^^ + 1] + C/ 



£2 
C 



In 



^r^^J[{B - A)^^ - 1][{B + A)i;^ + 1] - C/V^ 
-^^[{B - A)i>^l; - 1][{B + A)^^ + 1] + C/V^ 



B-A 
1+Aipip 



^ [{B - A)i,i^ - 1] [{B + A)i,i; + 1] - C/ V=3A 



(53) 
(54) 

(55) 



3 Identification of time in the classical model 

In the definition of an intrinsic time the Hamiltonian constraint plays a crucial role. In any constrained 

system we can impose the constraints in different steps. In classical mechanics, for example, we may 
first solve the equations of constraint to reduce the degrees of freedom of the system and obtain 
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a minimal number of dynamical variables. On the other hand, we may multiply the constraint by 
a variable parameter and add it to the Lagrangian. This Lagrange multiplier plays the role of an 
additional dynamical variable and the equations of motion consist of those obtained from variation 
of the Lagrangian with respect to the dynamical variables plus the equation of constraint. Solving 
this system of equations of motion leads to either the time evolution of the system or the forces 
of constraint. Also, when quantizing the system, we may impose the constraint before or after the 
quantization has been done. Now, if our system is the entire universe, e.g. in the case of quantum 
cosmology, these procedures result in different approaches to the problem of time rcparameterization. 
Here our goal is to reparameterize time in the classical model in a manner described below. We first 
solve the equation of constraint to obtain a set of genuine canonical variables with which to construct 
the Hamiltonian. Equations of motion are then obtained from this Hamiltonian and describe the 
evolution of the system with respect to this intrinsic time. 

To start, let us rewrite Lagrangian (38) in the following form 



L = upu + vpv + wpw + p^ip + ippj, + 



1 

12 



e-^-{pi-pi-pi)+e'-{A-V) 



3u/ 



(56) 



The role of the lapse function as a Lagrange multiplier is now clear in the above form of the Lagrangian. 
The term in the square brackets is indeed the Hamiltonian constraint. The procedure one should 
follow is the same as that described in [17]. First we make a choice of time in terms of one of 
the dynamical variables in the model and then solve the constraint equation (41) for its conjugate 
momentum and substitute the result in the Lagrangian. This process leads to a reduced Hamiltonian. 
A natural choice for time in terms of the intrinsic geometry would be the use of the scale factors. We 
may also take their conjugate momenta or a scalar character of the matter field. These latter choices 
are often referred to as extrinsic time. 

Let us now examine the above procedure in the following cases. First, suppose we choose t = u 
as an intrinsic time. Solving the constraint equation (41) for pu and substituting the result in (56) 
yields the square-root Hamiltonian 



H 



(2 2 



+ 12Ae' 



,6t 



(57) 



Here, we have assumed that the spinor field is free with mass m, i.e. V{'ip,ip) = mipil) and also used 
the relation ■0^' = 4e~^"pi^p^. Since this Hamiltonian is independent of v, w, i/j and jp, these variables 

are cyclic and their conjugate momenta are constant, that is Pv = ci, Pw = C2, p^ = C3 and p^ = C4. 
The functional form of all dynamical variables v, w, ^ and ^jJ resulting from (57) becomes 



1 



In 



VC±12|A|e«* + \/C 



(58) 



where = c{+(^ — 48mc3C4, with the upper and lower signs denoting positive and negative cosmo- 
logical constants respectively. A problem related to Hamiltonian (57) is that it is a time-dependent 
function. Such Hamiltonians describe a system which exchanges energy with the surrounding en- 
vironment. However, in the case of cosmology where the system under consideration is the whole 
universe, a surrounding environment does not have any meaningful interpretation. Therefore, such 
a Hamiltonian and the corresponding time parameter do not seem to be suitable unless the cosmo- 
logical constant is zero. It can be easily shown that in this case all dynamical variables have a linear 
behavior with time. Another feasible alternative for intrinsic time would be to take t = v. Following 
the same procedure described above we find 



H=(pI-pI- 12Ae^" + A^mp^p^) 



1/2 



(59) 



which does not suffer from being dependent on time. This Hamiltonian shows that the variables 
and -0 are cyclic and their conjugate momenta are constants, that is Pw = ci, Pip = C2 and p^ = C3. 
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The corresponding cosmology is then obtained from the solution of the equations of motion given by 
ii = etc. The result is 

n{t) = In (^e-'-^-^-^?j '\ (60) 
wit) = %t + wo, (61) 



C 

m = '-^t + ^o, (62) 



where uq, wq and V'o are integrating constants and = 48mc2C3 — . The above expressions show 
that this choice of the intrinsic time is appropriate for a universe with positive or zero cosmological 
constant. In the case when A = all canonical variables become cyclic and their time evolution 
becomes linear with time, as in equations (61) and (62). Choosing w as time gives rise to the same 
results. 

As we mentioned before, there are also time variables which are conjugate to the intrinsic time 
known as the extrinsic time. For example, consider the case where t = Pu- Solving the constraint 
equation for u leads us to the following Hamiltonian 



H 



{t^ -pI-Pw + '^^mp^p^ ' . (63) 



Vl2A 



This Hamiltonian is again time-dependent and describes a universe with positive cosmological con- 
stant. Since the variables v, w, tp and -0 are all cyclic their corresponding momenta are constants, 
give by Pv = ci-, Pw = C2, Pi) = C3 and p^ = C4. It can be shown that the evolution of all dynamical 
variables has the functional form 



(64) 



where = 48mc3C4 — cf — c^. Since the constraint equation (41) is independent of v and w the 
choices p„ and do not make good as time parameters in our model. In addition to the types of time 
described above, there is yet another choice of the extrinsic time which corresponds to the matter 
fields in the theory. In the case of a scalar field, the natural choice for time is the scalar field (j) itself, 
as has been done in [17]. In our model the matter field is a 4-spinor ^/^ with complex components 
and the simplest scalar constructed from it is ■j/''^. Indeed, it is easy to see that this quantity is 
proportional to the energy density of the spinor field, {pd = mtpci'il^cl)-, and is the only observable that 
enters in the classical Einstein equations. Thus, as the final time identification we take t = tpip and 
are led to the following Hamiltonian 

H = (l2Ae6« -pI+pI+ pI) ^'^ ■ (65) 



m 



The classical cosmology resulting from this Hamiltonian can be read from the equations of motion 
ii = etc., with the result 



«(*) = ln||^^i + e-^"°] , (66) 

-1 

, Pv = Cl, Pw = C2, (68) 
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where ci, C2, vq, wq and uq are integrating constants and = + It is clear from the above 
equations that this type of time is suitable for a universe with positive or zero cosmological constant. 
Again in the case of a zero cosmological constant the dynamical variables have a linear behavior with 
time. 

Our classical investigation on the problem of time is now complete. In the next section we shall 
pursue this problem within the context of a quantum cosmological model. 



4 Quantization of the model 



The usual approach to canonical quantization of a cosmological model is the Wheeler-DeWitt ap- 
proach where one uses the Dirac method to quantize the degrees of freedom of the system. The role of 
constraints in their operator form is to annihilate the wave function of the universe. This procedure 
leads one to the basic equation of quantum cosmology, the so called WD equation. This approach 
to quantum cosmology has its own problems and is not the subject of study in this paper, see [12] 
and [16]. However, as was done in the previous section, one may solve the constraint before using it 
in the theory, in particular before quantizing the system. If we do so, we are led to the Schrodinger 
equation 

(9* 

H^ = i-g^, (69) 

where H is the operator form of the reduced Hamiltonian, i.e. the operator form of one of the 
Hamiltonians given in (57), (59), (63) and (65). There are, however, some problems related to these 
Hamiltonians at the quantum level. Firstly, they are all given as square-roots and thus by means 
of the spectral theorem are assumed to be positive definite and Hermitian operators. Secondly, the 
Schrodinger equation (69) is in general a time-dependent equation because of the dependence of H 
on t. Thus one should not conclude that this equation is equivalent to a WD or a second order 
Klein-Gordon type equation like 

H'^ = (70) 

to remove the square-root form of the Hamiltonian. Indeed, by acting with H on both sides of 
equation (69) one obtains 



which, of course, has different solutions from those given by equation (70), unless the Hamiltonian 

is timc-indcpcndcnt. For a careful description of the details of this issue see [14, 17]. However, as 
is well known from the elementary quantum mechanics, if the Hamiltonian is time-independent, the 
solutions of the Schrodinger equation can be written as 



^{x,t) = exp [-iH{t - to)] '^{x,to). 



(72) 



In the case of a time-dependent Hamiltonian the above formula should be modified by the Dyson 
series 



*(x, t) 



°° ft fti ft„-i 



In practice, interesting cases happen when 

[H{t'),H{t")]=Q, 

in which case the Dyson series becomes 



^(f,to). (73) 



(74) 



exp \-i H{t')dt' 



*(x, to). 



(75) 
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Another feature of relation (74) is that the Hamiltonian has the same cigcnstates at all times, in 
other words if ^'eq(x) is the eigenstate of H(to) at some initial time to : H{to)"i/ Eo{x) = Eo'^Eoix), 
then '^Eo{x) is also the eigenstate of H{t) with another eigenvalue E{t) at time t 

HmEo{x)=EmEoix). (76) 
In conclusion we are led to the following relationship for the time evolution of the wave function 



exp(^-i^ E{t')dt' 



^E^to)- (77) 



To find \I'Eg(x, io) we can use the eigenvalue equation //(to)^_Bo(^) = -E-O^Eol^) ™ the form 

H\to)^Eoix) = Ei^Eo{x). (78) 

This is an allowed equation since H is evaluated at a particular fixed time to Siud thus the square-root 
problem is now resolved. In view of dealing with the Klein-Gordon or WD equation in our Schrodinger 
approach, we do not need to be concerned with the inner product of states since the determination of 
the wave functions in the Schrodinger equation and the construction of the Hilbert space proceed in 
the usual sense of quantum mechanics. The final remark about the rcdiiced Hamiltonian in the above 
procedure is the factor-ordering problem when one embarks on constructing a quantum mechanical 
operator equation. In the class of Hamiltonians represented by (57) and (63), this problem is not too 
important since there is no presence of a canonical variable with its conjugate momentum. However, 
consider for example, Hamiltonian (59) in which both u and its conjugate momentum p„ are present 
and do not commute. In dealing with such Hamiltonians at the quantum level care should be taken 
when one tries to make a change of variable. This is so because such an operation may change the form 
of the Hamiltonian where a product of a variable with its conjugate momentum has now appeared. 
This is an indication that in quantizing the system, the ordering problem becomes important. As an 
example, consider a change of variable = e^" applied to Hamiltonian (59) , changing it to 

/9 \^/2 
H=i -r^Pr -pI + ^^mp^p^ - 12ArM . (79) 



Under the same change of variable Hamiltonian (65) becomes 



H 



1/2 

(80) 



It is now clear that in replacing r and p^ with their corresponding operators, the ordering considera- 
tions should be taken into account. We return to this issue again in the next section. 



5 Time identification in the quantized model 



Let us now use the above theory in the problem at hand and investigate time rcparamcterization 
introduced in the last section in the context of the quantized model. In the case when time is chosen 
as t = u, the Hamiltonian is given by (57) satisfying the relation (74) for all times. To find the 
wave function ^'(■u, lu, t), we must first solve equation (78) for this Hamiltonian. With the usual 
replacement Pv ^ and similarly for p^, p^ and p^ this equation becomes 



92 



- — - + 12Ae^*o + 48m 



92 



^Eoiv,w,tpip) = Eo^Eoiv,w,tp'4;). 



(81) 



The solutions of the above differential equation are separable and may be written in the form 
'^{v,w,ijjilj) = V {v)W (w) f (iptp) leading to 



1 £v 



a 



1 £w 

W dw^ 



(82) 
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1 d'f 



7 



(83) 



where a, j3 and 7 are separation constants satisfying + (5"^ = 48m7^ + 12Ae^*° — Eq. Equations 
(82) have simple solutions in the form of exponential functions e~l"l^ and e'~l^l"'. The exponents are 
chosen so that ^{v,w — +00) = 0. To find the solutions of equation (83) we use the ansatz 



(84) 



n=0 



After a little algebra we find 



(n!) 



Co and 



00 2n 

/W) = E7^W)" 



=0 (^0 



(85) 



It is easy to check that this series converges for all value of iptp. We may thus write the solution of 
equation (81) as 

-\a\v- 



00 2n 
■\w / 



E 



(86) 



n=0 ^ ' 

According to equation (76) the above eigenfunctions should also be the eigenfunctions of H{t) given 
by (57) with eigenvalues E{t) such that 



E{t) 



12A (, 



)+El 



1/2 



(87) 



The time evolution of the wave function is then given by (77) with the result 



exp 



to 



E{t')dt' 



jexp (-i £ [l2A(e6*' - e^*") + E^] dt'^ | ^Eo{v, w, 



= jexp 



72A 



2 [l2A(e'^* -e'^*o) + £;^ 



1/2 



-2Eo + {E^-12Ae^^°): 



In 



{ [12A(e6* - e6*o) + E^]^^'' - {El - 12Ae6*o)V2| ^^Eq + {El - 12Ae6*o)V2| 



{ [l2A(e6* - e6*o) + El] ^''^ + (^2 _ i2Ae6to)i/2| [Eq _ (^2 _ i2Ae6*o)i/2} 

'J'£;o(i;,W,'0^/;). 

A comment about the Hamiltonian (57) is that its eigenenergies have a continuous spectrum in the 
range [\/l2Ae^*, +00) which for, A > 0, is positive definite as required. This means that our choice of 
the time parameter is suitable for a universe with positive or zero cosmological constant. However, 
in the case of A = 0, the time evolution of the wave function is simply given by (72). Note that in 
the classical model in this gauge both positive and negative cosmological constants were valid. 

Let us now deal with the choice t = u as time. In this case the Hamiltonian is given by (59). 
However, to deal with the ordering problem we work with Hamiltonian (79). To guarantee Hermiticity, 
the operator form corresponding to this Hamiltonian should be written as 



pIj + 48mp^p^ - 12Ar*^ 



1/2 



(89) 



where the parameter p denotes the ambiguity in the ordering of factors r and Pr in the first term of 
(79). Taking p = 0, equation (78) reads 



9d_ 2 a_ 

4 dr dr dw^ 



48m 



92 



12A/ ^Eo{r,w,'ilji;) = eI^ Eo{r,w,ipi^). 



(90) 
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We again write the solutions in the form "^{rjWjipip) = R{r)W (w) f (ipip) and find 



1 d'^W 2 



(91) 



W dw^ 

P\ (92) 



1 1 «2 



/ dtjjdtfj 48m 

where a and P are separating constants and 7^ = — The solutions of equations (91) and (92) 
are known as 

W{w) = e-M^ and /(^^) = E ^^^f^^^'" W)" 

n=0 

Also, the solutions of (93) for a positive cosmological constant can be written in terms of Bessel 
functions as 

i?(r)=r-V2j±,QV3Ar2), (94) 
where i/^ = 1/16 — {Eq — 7^)/9. The complete solution of equation (90) now reads 



where to avoid diverging solutions at r = we have removed the function J^^,. As is clear from the 
above equation the wave function satisfies ^{w +00) = and ^{r +00) = 0. In the limit r — ^ 
the wave function ^'(r, it;, ^'(/') behaves as r'^'^~^/'^ and thus to have regular solutions near r = we 
must have 2^ — 1/2 > 0. This condition restricts the eigenenergies to the interval [0, 7] and the initial 
condition to *(r = 0) =0 for the wave function as suggested in [1]. Contrary to the classical case 
where such a choice for the time parameter was appropriate only for a positive cosmological constant, 
A can be negative as well as positive here. In the case of a negative cosmological constant, the Bessel 
function in equation (94) must be replaced with the modified Bessel functions Ki,{x) and Iv{x). To 
satisfy ^{r ^ 00) = we restrict ourselves to functions Ki,{x). Thus for a negative cosmological 
constant 



, 2n 



/2 I \ °° f/3/\/48mV 

^E,{r,wM) = e-\^\^r-^I^KJ-p^\rAY, ^ / {i^^l^T . (96) 

Again to avoid singularity at r = the order of the function Ki,{x) should be pure imaginary; i/^ < 
[32], which results in the interval {\/^'^ + 9/16, +00) for the allowed eigenenergies. One should note 
that it is impossible for the wave function to satisfy the condition ^'(r = 0) = in this case. 

For a universe with zero cosmological constant the solutions of equation (93) can be written as 

R{r) r''+ .r"- , (97) 

where 



'^ = 2 



{-\±^i-\Q{Ei-^^)/^y 



Thus for the energies in the interval < £^0 < 7 we can have regular solutions r'^+ satisfying ^{r = 
0) = 0, while for \/^'^ + 9/16 < Eq < +00 we are led to oscillatory solutions. Since Hamiltonian (89) 
is time-independent the time evolution of ^'(x) is given by (72) with the result 

*(r,'w;,'0V',t) = e-'^°^*-*°^^Eo{r,w,tj;^). (98) 
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The other choice for time studied in the last section is t = which may be interpreted as the 
conjugate to the choice t = u. It can easily be shown that in this case the eigenstates of Hamiltonian 
(63) are obtained by interchanging the roles of to and Eq in equation (88) [17]. 

Our final discussion about the problem of time is to take a scalar character of matter as the time 
parameter. As is well known the simplest scalar which can be constructed out of a spinor field is -ipip. 
This suggests the choice t = tjjtp. In this case we must solve equation (78) with Hamiltonian (65) or 
(80) to deal with the factor ordering. With the same factor ordering parameter as we have taken in 
equation (89) we write 



dr dr 



+ 



9 \ dv"^ dw"^ 



+ 



16A 



*B(,(r,i;,u;) 



64 



mEl-^Eo{r,v,w), 



(99) 



with solution 



^ ^-3/2{\a\v+\i 



where u"^ = 1/16 + (64m£^g + 3-f'^)/12 and 7^ 



1-)^-l/27 



(100) 



a 



+ /J^. The same analysis done for equation 



(95) suggests that if A > for all energies in the interval [0, +00), we have solutions satisfying the 
condition ^'(r = 0) =0. In the case of a negative cosmological we get solutions that are neither 
regular as r ^ nor in the limit r — >■ 00. Thus, as in the classical case this choice of time is only 
suitable for a universe with A > 0. 



6 Conclusions 

In this paper we have studied the time reparamctcrization problem in a Bianchi type I cosmology 
with a Dirac spinor field as the matter source at both the classical and quantum levels. This problem 
arises from the fact that Einstein field equations in classical general relativity are under-determined, 
resulting in the requirement of imposing a gauge condition before a solution can be found. In the 
ADM formalism of general relativity these gauge conditions are those that are commonly imposed on 
the lapse function or shift vectors. 

With line element (1) we have taken the shift vectors A^* = so that the only gauge freedom in 
our model relates to different choices of lapse function N{t), giving rise to different choices of time 
parameters in the model. At the quantum level however, the gauge freedom appears in the form of the 
Hamiltonian and momentum constraints. Again, because of the form of our metric the momentum 
constraint is automatically satisfied and one should only deal with the Hamiltonian constraint. In 
order to fix the gauge, we chose a time parameter defined in terms of the 3-geometry (intrinsic time) 
or its conjugate momentum. We also explored the possibility of the matter field playing such a role 
(extrinsic time). The Hamiltonian constraint for the conjugate variables was then solved and the 
solutions, considered as time, were used to find a reduced Hamiltonian and consequently the classical 
field equations. Of course, these equations do not suffer from the under-determinacy problem. In 
this sense we have taken three intrinsic times as t = u,v,w where u, v and w are related to the scale 
factors through equations (37). The corresponding cosmologies in terms of these time parameters 
are given by the relations (58) and (60)- (62). These solutions are valid for an arbitrary cosmological 
constant when t = u and a positive or zero cosmological constant in the case oft = vovw respectively. 
Also, the extrinsic times, t = Pu and t = iptp, give rise to the classical cosmologies (75) and (66)-(68), 
describing a universe with a positive, negative or zero cosmological constant respectively. We have 
also shown that py and Pw are not suitable for playing the role of time. To quantize the model we 
have followed the procedure introduced in [17]. After fixing the time we used the operator form of the 
reduced Hamiltonian to quantize the system in the Schrodinger picture. This procedure has led us to 
a quantum cosmology with time dependent wave function (88) in the case t = u which is appropriate 
for A > 0, or (95)-(98) when t = v which is suitable for a cosmological constant with arbitrary sign. 
Finally, in the case where -t/jip was considered as the time variable we have obtained the wave function 
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(100) for a positive or zero cosmological constant. 
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